-1- 



ACTIVE INTERFEROMETRIC SIGNAL ANALYSIS IN SOFTWARE 

Field of the Invention 

The invention generally relates to active signal analysis systems, and in particular, 
5 interferometric signal analysis systems. 

Related Applications 

This application is a continuation-in-part of U.S. Pat. App. No. 09/913,132, filed 
Aug. 9, 2001, entitled "Method and System for Signal Detection in Arrayed 

10 Instrumentation Based on Quantum Resonance", which is a national stage application of 
PCT Application No. USOO/04076 filed on Feb. 17, 2000, which claims priority from U.S. 
Pat. App. No. 09/253,789, now U.S. Pat. No. 6,136,541, filed on Feb. 22, 1999, entitled 
"Method and Apparatus for Analyzing Hybridized Biochip Patterns using Resonance 
Interactions Employing Quantum Expressor Functions, all of which are hereby fully 

1 5 incorporated by reference". Also incorporated herein by reference in their entirety are the 
following patents based on applications filed contemporaneously with U.S. Pat. App. No. 
09/253,789: U.S. Pat. No. 6,142,681, entitled: "Method and Apparatus for Interpreting 
Hybridized Bioelectronic DNA Microarray Patterns Using Self-scaling Convergent 
Reverberant Dynamics" and U.S. Pat. No. 6,245,511, entitled "Method and Apparatus for 

20 Exponentially Convergent Therapy Effectiveness Monitoring Using DNA Microarray 
Based Viral Load Measurements". Also incorporated herein by reference in their entirety 
are the following provisional patent applications: U.S. Pat. App, No. 60/395,074, filed on 
Jul. 9, 2002, entitled "Active Signal Processing Technique Using Resonance 
Interferometry". U.S. Pat. App. No. 60/399,227, filed on Jul. 29, 2002, entitled "High- 

25 Throughput Quantitation of Signals from Arrayed Platforms Using Resonance 
Interferometry" U.S. Pat. App. No. 60/408,057, filed on Sep. 3, 2002, entitled "Technique 
for Analysis and Quantitation of Arrayed Information Structures". U.S. Pat App. No. 
60/408,844, filed on Sep. 5, 2002, entitled "System and Method for Analysis and 
Quantitation of Arrayed Platforms". U.S. Pat. App. No. 60/440,253, filed on Jan. 15, 2003, 
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entitled "Technique for Detection and Quantitation of Arrayed Platform Signals", 

Background of the Invention 

L Overview of Passive and Active Interferometric Signal Analysis 
5 Interferometric analysis has proven useful for performing signal analysis, 

particularly for detecting events of interest within signal patterns, such as within a two- 
dimensional image. One example of interferometric analysis involves irradiating a moving 
object, such as an aircraft, with radar pulses, that interact with the object and interfere with 
one another and which are then reflected back to a sensing device such as an antenna. 

10 Signals received by the antenna are compared against signals originally transmitted in an 
effort to detect a unique spectral signature of the object so as to identify the object. This 
form of interferometric analysis is referred to as "active analysis", because signals used to 
irradiate the event of interest interact with the object and with other signals to actively 
modify the resulting interference pattern. In other words, the system being analyzed is 

1 5 actively modified. In the case of radar, electromagnetic fields surrounding the event of 
interest are actively modified. 

An alternative type of signal analysis is "passive analysis" wherein signals that 
have already been detected are processed in an effort to identify events of interest therein. 
Examples include filtering images in an attempt to enhance features therein. These 

20 techniques are deemed to be passive because no noise or other signals are added into the 
system being analyzed. More specifically, passive interferometric analysis tj^ically 
involves the manipulation of observed intensity patterns or phase information without the 
addition of noise or other signals. Often, passive analysis systems exploit parametric 
methods such as matched filter techniques, maximum a-posteriori (MAP) techniques, 

25 maximum likelihood estimator (MLE) techniques, singular value decomposition (SVD), 
classical Fourier analysis, parametric distributional clustering (PCA), or median filter 
variants. Non-parametric methods may be employed as well. In addition, quantitation and 
clustering techniques may be employed such as "soft-computing" and its variants which 
include neural networks, fuzzy logic systems, and Bayesian inferencing systems. 
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Although a wider range of passive analysis techniques have been exploited, each requires 
that the observed signal have an intensity greater than the intensity of the background 
signal in which it appears. In addition, the analysis systern must properly estimate the 
background signal and to properly estimate any variability in the background signal In 
5 some cases, spread spectrum (S-S) amplitude correlation is provided with background 
subtraction, and scaling and normaUzation is employed to compensate for detector 
nonlinearities. Hence, there are various limitations to passive analysis, particularly the 
inability to detect signals having an intensity less than or near to the intensity level of the 
background, 

10 Insofar as the analysis of hybridized biological microarray patterns is concerned, 

heretofore, passive analysis has been the only form of analysis performed. Briefly, a signal 
pattern output from a biological microarray is normahzed to compensate for dye, scanner, 
or array effects . Then an estimate is made as to the background signals at feature, 
transcript or chip levels. The background is then subtracted from the foreground and a 

15 ratio of foreground minus background is calculated to qUantitate a differential expression 
representative of patterns of interest. Often, this calculation is performed using channel 
ratios. In any case, examples of passive analysis of hybridized digitized image patterns 
include PCA, SVD, MLE, hidden Markov Model (HMM), and various statistical 
approaches. As the technique is passive, it is subject to the Hmitations summarized above, 

20 in particular, only hybridized transcripts with genes above the estimated background can 
be detected. In addition, high call variability in targets with low expression levels arises 
and high CV in fold change computations occurs. Skewed results can arise when dealing 
with platforms or protocols having a limited dynamic range. Typically, passive 
hybridization analysis techniques require vast amounts of data and exemplars to allow 

25 refinernent of statistics driving the call analysis. Heretofore, efforts to improve the 
accuracy and reliability of passive techniques for analyzing hybridized digitized images 
involve wet signal amplification and biological enhancements. Computational 
enhancements of the passive analysis have been fairly minimal. 

Thus, heretofore, computer systems have been largely limited to performing 
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passive analysis whereas active analysis is largely performed only by hardware systems. 

2. Hardware-Based Active Signal Analysis 

With hardware-based active interferometric signal analysis, external excitation is 
provided to the system being analyzed while the system is being analyzed to provide an 
5 interaction that generates additional information or provides some change in a qualitative 
property of the system. Typically, the external excitation involves some form of noise 
injection, which results in a nonlinear coupling between the noise and events of interest 
within the system being analyzed, which enhances those events. As noted, one example of 
active signal processing (ASP) is radar wherein the signal being injected into the system 

10 being analyzed is the radar signal Another example is the superconducting quantum 
interferometric device (SQUID), which utilizes the quantum effects of superconductivity to 
detect magnetic fields. Other examples include imagery intelligence (IMINT), signals 
intelligence (SIGINT), electronic intelligence (ELINT) devices. In addition, active signal 
processing has been applied to femto-lasers. While these hardware-based systems provide 

1 5 enhanced signals in some circumstances, such systems are expensive to build and operate 
and can only physically manipulate or generate very specific types of signals. 

3. Stochastic Resonance Techniques 

Another exemplary active system or technique that has been implemented in 
20 software is stochastic resonance wherein noise is applied to an image or a system being 
analyzed to improve at output signal-to-noise ratio (SNR). With stochastic resonance, a 
detector is modeled in terms of a bistable nonlinear dynamic system, which is enhanced by 
applying additive random noise element and a . periodic sinusoidal forcing function. 
Stochastic resonance occurs when the SNR passes through a maximum as the noise level is 
25 increased. By properly selecting the amount of noise being applied, events of interest can 
be enhanced so that signals below background can be detected. 

Traditionally, stochastic resonance is employed using a classical noise source such 
as a standard Gaussian noise pattern. Alternatively, however, stochastic resonance may 
employed using some form of quantum noise. In this regard, it has been shown that 
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irradiation of a bistable system with a weak resonance laser field induces long-lived 
coherent tunneling motion whose amplitude exhibits a quantum stochastic resonance 
(QSR) as a function of the two-level-system relaxation rate and the field strength (See 
Makarov DE. Makri N. "Stochastic resonance and nonlinear response in double-quantum- 
5 well structures" Physical Review B, 52(4):R2257-R2260, 1995 Jul 15). Essentially, QSR 
is an enhancement of the response of a driven quantum mechanical system by quantum 
noise. 

In any case, heretofore, traditional stochastic resonance and QSR have only been 
applied to point dynamical systems, i.e. to systems providing only a signal data point 
10 varying in time. One other technique for enhancing stochastic resonance is so-called array- 
enhanced stochastic resonance (AESR), which extends stochastic resonance to 2-D. 
dynamical systems, though still limited to static systems. Although an improvement, 
AESR is typically limited to the analysis of systems with relatively few degrees of 
freedom, 

15 Thus, most conventional active signal analysis systems are not implemented in 

software and those that have been implemented in software typically have significant 
limitations. Accordingly, it would be desirable provide new or improved techniques for 
implementing active signal processing in software and it is to these ends that aspects of 
invention are drawn. By implementing active signal analysis in software, all of the 

20 advantages associated with software systems over hardware systems can be exploited, 
including lower cost, increased flexibility, increased robustness, portability, etc. 

More specifically, it is desirable to provide an active signal processing system 
implemented in software capable of detecting events of interest in spatial, temporal or 
spatio-temporal systems within a combinatorial state space wherein the state space is 

25 discrete and continuous and wherein background noise is nonstationary. Ideally, the 
system would be capable of analyzing signals received from a wide variety of arrayed 
platform detector configurations including linear, arrayed, spectral, and time-varying and 
can handle TV-attribute systems. Moreover, ideally, the improved system would be capable 
of analyzing signals with a signal to background ratio of 10"^. Insofar as the analysis of 
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biological samples applied to a microarray platform is concerned, the system or technique 
would ideally be capable of detecting and quantitating specific hybridization (e.g., genomic 
events) and specific binding events (e.g., proteomic events). Insofar as genomics is 
concerned, the system should be capable of processing data from a high-density gene 
5 expression microarrays (i.e. array is providing 10^ to 10^ spots per array) and which 
represents discrete, spatially indexed systems having high probe densities (i.e. from 10^ to 
10^ probes per spot). 

Summary of the Invention 

10 In accordance with the invention, techniques are provided for performing active 

interferometric signal analysis in software. Depending upon the implementation, the 
techniques provide for both detection and quantitation analysis by exploiting constructive 
or destructive interferometric analysis (or a combination of constructive and destructive 
interferometric analysis) usiiig reverberant convergence to detect resonance events. The 

1 5 techniques achieve software erhulation of wave-particle interactions and wave-wave 
interactions and can operate in either the frequency domain or the phase domain. The 
techniques may be used for analyzing static spatial systems, static data from arrayed 
measurement platforms, dynamical systems, spatio-temporal systems or plasma systems. 

Briefly, the techniques exploit expressor functions designed to reject any 

20 interfering noise or background clutter from possible events of interest within arrayed data 
to be analyzed. More specifically, the expressor functions are designed so as to extract 
spectral invariants of possible events of interest associated with the array device used to 
detect the arrayed data. Various aspects of the invention are directed to the generation of 
the expressor fiinctions. Typically, the expressor functions are designed to incorporate 

25 some form of non-Gaussian noise, such as quantum noise, pseudorandom noise, 
polycyclostationary noise, cyclostationary noise, stationary noise, non-stationary noise or a 
systemic bias. Preferably, prior to the application of the expressor function to arrayed data, 
the array is pre-conditioned so as to convert the arrayed data to a spectral domain in which 
spectral hamionics parameterize events of interest to a pre-deterrnined dynamical system. 
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Preconditioning may be performed by applying a preconditioning function incorporating, 
e.g., a 1-D Fourier function, a 2-D Fourier fiinction, an N-D Fourier function, a time 
division multiplexing (TDM) fiinction, a wavelength division multiplexing (WDM) 
function, a frequency division multiplexing (FDM) function, a radial bias function, a 
5 wavelet kernel function, a fractal function, or a soli ton function. The arrayed data itself 
may be the form of a spatial 1-D array, a spatial 2-D array, a spatial point emitter array, a 
temporal point emitter array, a spatio-temporal point emitter array, a spectral point emitted 
array or a virtual array constructed by combining spatial separate point emitters. 

In accordance with an exemplary method, a technique is provided for actively 

10 analyzing a signal pattern representative of arrayed data to identify events of interest 
therein. A signal pattern representative of arrayed data is input and resonance patterns are 
generated based on interference between synthetic noise and the signal pattern. Next, 
resonances within the resonance patterns associated with events of interest are detected. 
Preferably, the signal pattern is preconditioned prior to the step of applying synthetic noise 

15 to the signal pattern, 

In a first exemplary system embodiment, also referred to herein as an "open loop 
interferometric system", a system is provided for analyzing an arrayed signal pattern 
generated by an arrayed platform device to identify events of interest within the signal 
pattern. The system includes an expressor function input unit for inputting previously 

20 generated expressor functions capable of extracting spatial, spatio-temporal or spectral 
invariants of events of interest associated with the particular arrayed platform device being 
used. The system also includes a preconditioning unit for preconditioning the arrayed 
signal pattern. An active interferometric coupler operates to convolve the preconditioned 
arrayed signal pattern and the expressor functions so as to interferometrically enhance 

25 portions of the preconditioned pattern associated with events of interest, if any, present 
within preconditioned pattern. The system also includes a resonance marker detector for 
identifying the occurrence of events of interest within the convolved signal pattern. 

In a second exemplary system embodiment, also referred to herein as a "goal- 
directed interferometric systern", the system for analyzing an arrayed signal pattern 
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includes an expressor function input unit for inputting previously-generated expressor 
functions and an adaptive interferometric coupler for preconditioning the arrayed signal 
pattern so as to convert the arrayed signal pattern to a spectral domain while 
simultaneously convolving the signal pattern and the expressor fiinctions so as to 
5 interferometrically enhance portions of the signal pattern to identify events of interest, if 
any, within the signal pattern. 

In a third exemplary system embodiment, also referred to herein as a "self- 
organizing interferometric system", the system includes a preconditioner for 
preconditioning the arrayed signal pattern and an expressor function adaptation unit for 
10 generating preconditioned expressor fiinctions based on input canonical expressor 
functions (i.e. general expressor function not associated with the particular platform being 
used) and based on a preconditioned signal pattern while interferometrically enhancing 
portions the preconditioned signal pattern so as to identify events of interest, if any, within 
the signal pattern. 

15 In a fourth exemplary system embodiment, also referred to herein as an "iterative 

interferometric system", the system includes an expressor function input unit for inputting 
previously-generated expressor functions derived for the platform being used, a 
preconditioner for preconditioning the arrayed signal pattern and a iterative convolution 
coupler. The system also includes an adaptive controller unit for controlling the coupler to 

20 iteratively and selectively convolve the expressor functions to a preconditioned signal 
pattern until a predetermined degree of convergence is achieved so as to identify events of 
interest within the enhanced signal pattern. 

A variety of techniques are set forth for performing active signal analysis in 
software. Method, system and computer product embodiments are provided. 

25 

Brief Description of the Drawings 

The features, objects, and advantages of the present invention will become more 
apparent from the detailed description set forth below when taken in conjunction with the 
drawings in which like reference characters identify correspondingly throughout and 

r 
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wherein: 

FIG. 1 is a block diagram providing an overview of system components of the first 
exemplary embodiment of the invention (the open loop interferometric system), which 
provides separate expressor function generation unit, preconditioner, active interferometric 
5 coupler and resonant marker detector components; 

FIG. 2 is a flowchart providing an overview of method steps for the first exemplary 
embodiment of the invention; 

FIG. 3 is a block diagram providing an overview of system components of the 
second exemplary embodiment of the invention (the goal-directed interferometric system), 
1 0 which provides an adaptive interferometric coupler; 

FIG. 4 is a flowchart providing an overview of method steps for the second 
exemplary embodiment of the invention; 

FIG, 5 is a block diagram providing an overview of system components of the third 
exemplary embodiment of the invention (the self-organizing interferometric system), 
1 5 which provides a expressor fiinction adaptation unit; 

FIG. 6 is a flowchart providing an overview of method steps for the third 
exemplary embodiment of the invention; 

FIG. 7 is a block diagram providing an overview of system components of the 
fourth exemplary embodiment of the invention (the iterative interferometric system), 
20 which provides a adaptive controller for iteratively controlling a separate iterative 
interferometric coupler; 

FIG. 8 is a flowchart providing an overview of method steps for the fourth 
exemplary embodiment of the invention; 

FIG. 9 is a graph illustrating conibinations of arrayed data signal patterns that may 
25 be processed by any of the embodiments of the invention; and 

FIG. 10 is a graph illustrating combinations of preconditioning fiinctions and 
expressor functions that may be used in connection with any of the embodirnents of the 
invention; 

FIG. 11 is a graph illustrating an exemplary platform array detector and active 
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interferometric analysis system that may exploit the invention; 

FIG. 12 is a graph illustrating various platform array detectors that may be used in 
connection with any of the embodiments of the invention; 

FIG, 13 is a graph illustrating the structure of an exemplary microarray platform 
that may be used to generate arrayed data for use by any of the embodirnents of the 
invention; and 

FIG- 14 is a flowchart summarizing the overall method of the invention, at least in 
accordance with an exemplary embodiment thereof 
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Detailed Description of Exemplary Embodiments 

With reference to the remaining figures, exemplary embodiments of the invention 
will now be described. 

5 A. OPEN LOOP INTERFEROM ETRIC SYSTEM 
L Overview 

FIG. 1 illustrates, at a high-level, exemplary system components for the first 
exemplary implementation of the invention. Components employed only during a design 
phase are shown on the left; components employed during the actual analysis of an input 

10 signal are shown on the right Briefly, during the design phase, expressor function 
generation unit 100 generates one or more designer expressor functions based upon the 
specific characteristics of a platform array detector 102 to be used to detect signal patterns 
to be processed and also based on calibration data derived from platform array 102. The 
calibration data is generated based on known events of interest. That is, input signals or 

15 physical samples containing known events of interest are applied to the detector to 
generate the calibration data. In one example, the platform array detector is a genomic 
biochip/microarray and the calibration data is derived from a biological sample containing 
known gene expressions. The characteristics of the platform array detector that are input 
by the expressor lunction generation unit include both the layout of the array as well as its 

20 coherent noise or background signal characteristics. Once the expressor fimctions have 
been generated for the particular platform array detector based on the calibration data, the 
design phase is complete and, thereafter, the expressor function generation unit is no 
longer required unless an alternative platform array detector is employed, 

Components employed during the actual operational phase of the system include 

25 the aforementioned platform array detector 102 that generates a detected signal pattern and 
a preconditioner unit 104 that precotiditions the detected signal pattern so as to convert the 
signal pattern to a spectral domain (for an example of a technique for transforming the 
signal pattern to a spectral domain, please see U.S. Pat. App. Ser. No. 10/430,664 entitled 
"Method and System for Characterizing Microarray Output Data" filed on May 5, 2003 
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which is hereby fully incorporated by reference) in which spectral harmonics parameterize 
events of interest to a predetermined dynamical system (for which the applicable 
characteristics are input into the preconditioner unit. The operational phase components 
also include an active interferometric coupler 106 that convolves the preconditioned signal 
5 pattern and the previously generated expressor fiinctions. The convolution is performed so 
as to interferometrically enhance portions of the preconditioned signal pattern associated 
with events of interest, if any, that are present within the preconditioned signal pattern. As 
explained more fully below, the convolution is performed via reverberant convergence so 
as to emulate active interferometric enhancement. Finally, the system a;lso includes a 

10 resonant marker detector 108 that processes the convolved signal pattern so as to identify 
particular events of interest, if any, appearing within the convolved signal pattern. 

With the exception of the platform array detector, which is a physical hardware 
component, all other components illustrated in FIG. 1 may be implemented in software, 
hardware, firmware or some combination thereof In particular, each component may be a 

15 computer product operative to perform the functions described. In one example, each 
component is a software module operating within a single generally programmable 
computer. In other cases, the software modules are implemented using separate 
microprocessors or application specific integrated chips (ASICs). In still other examples, 
some of the components are implanted in software whereas others are completely hard- 

20 wired. As can be appreciated, a wide range of possible implementations is consistent with 
the invention and no attempt is made herein to describe all possible configurations. 

Method steps associated with the open loop interferometric system of FIG, 1 will 
now briefly be summarized with reference to FIG. 2. Although these steps are 
advantageously performed by the system of FIG* 1, the steps can be performed by any 

25 other suitable system. Conversely, the system of FIG- 1, or portions thereof, may be 
employed to perform methods different from that which is illustrated in FIG* 2. The same 
is true for the various other figures of the patent application. In any case, at step 1 1 0, a 
signal pattern derived from a particular platform array detector is input and, at step 112, 
expressor functions designed to extract spectral invariants for events of interest detectable 
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by the particular platform array are input. Then, at step 114, the input signal pattern is 
preconditioned to convert the signal pattern to the spectral domain in which spectral 
harmonics parameterize events of interest to the aforementioned predetermined dynamical 
system. At step 116, the preconditioned signal pattern and the expressor functions are then 
5 convolved using reverberant convergence so as to emulate an active interferometric 
enhancement of portions of the signal pattern associated with events of interest thus 
yielding a convolved signal pattern (also referred to herein as an enhanced signal pattern. 
Finally, at step 1 18, the convolved signal pattern is examined to identify events of interest 
therein. 

1 0 It should be noted that a list of specific events of interest need not be input at step 

1 12. Rather, the technique operates to detect possible events of interest present in the input 
pattern based solely on the input pattern and the previously-generated expressor functions. 
Hence, the technique of FIG. 2 is not merely a diagnostic technique that seeks to 
determine whether particular pre-determined events of interest are present in the input 

15 signal pattern. This will be more readily apparent from the detailed examples which 
follow. 

2. Exemplary Implementation 

The mathematical equations performed by the various components of FIG. 1 or 
20 during the various method steps of FIG. 2 will now the set forth in detail In this example, 
quantum expressor functions are employed. However, as described below, other expressor 
functions can alternatively be u$ed, with the equations modified as needed dependent upon 
factors such as signal to background strength, extraction core length of the events of 
interest being analyzed, platform coherence, and desired sensitivity. 
25 The exemplary equations are particularly relevant for conditions where the 

background is 10 to 1000 times stronger than the signal. 

Initially, during the design phase, expressor function generation unit 100 of FIG. 1 

inputs platform array detector characteristics represented mathematically as: 
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where Q{N ,M\T,A) denotes a preconditioned extraction core corresponding to 
one or more events of interest. The indices A'^ (0:n), ]Vl(0:m), T(l:k), A(0:p) respectively 
refer to the spatial, temporal and spectral dimensions of the physical or virtual arrayed 
platform described in FIG 11, Details of preconditioning the extraction core corresponding 
to event of interest are given in U.S. Pat, App. Ser. No, 10/430,664 entitled "Method and 
System for Characterizing Microarray Output Data" filed on May 5, 2003. 

The Expressor Function Generation unit also inputs calibrated data (j exemplars) 
represented mathematically as: 



Based on the input calibrated data, the expressor function generation unit performs 
calculations provided below. 

In some embodiments. Quantum Expressor Functions (QEF's) are calculated based 

upon the signal pattern based upon one or more events of interest. The QEF is generated 

by the steps of: 

(a) selecting a canonical coupled Spin Boson basis system by first 



conditioning a Haniiltonian for that system; 



(b) 



calculating harmonic amplitudes |Pml for the Hamiltonian 



(c) 



generating an order function (OF); 



(d) 



measuring the entrainment states of the OF using exemplars with ground 



truth data; and 



(e) 



modulating the OF using calibrated exemplars to yield the QEF 
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FuUy characterized exemplars of the event of interest are used to develop a 
platform specific QEF/ Hamiltonians are utilized as they provide a mathematical basis for 
determining the dynamical behavior of a system at each coordinate and at each momentum. 
5 Hamiltonian systems are further characterized by kinetic and potential energy. A key 
property of Hamiltonian systems is that they are conservative and have no dissipation of 
energy during relaxation. The use of a Spin Boson Hamiltonian, which is a particular type 
of Hamiltonian systems, permits the exploitation of quantum stochastic resonance 
phenomenonology. Hamilton's Equations and Hamiltonian ftmctions pertain to well 

1 0 understood concepts in dynamical systems theory and dynamical equations of motion. (See 
"Theory and Problems of Lagrangian Dynamics", Ed. Dare A, Wells, Schaum*s Outline 
Series, McGraw-Hill Company, NY, 1967, pages 316-321; Elements of Hamiltonian 
Mechanics", Dter Haar, Pergamnon Press, 1960, pages 95-102; and "Statistical 
Mechanics", R.K. Pathria, Pergamon Press, Toronto, 1972, pages 6-7, 136-147). 

15 The QEF encodes quantum noise to drive order in an asymmetric bistable 

quantum-mechanical system. The appearance of a resonance requires an asymnietry in the 
energies of the two states. A rate equation constructed for transitioning between the two 
states of a bistable system, such that the dynamics can be characterized in terms of 
transition rates O H- and O-^ between two asymmetric quantum superposition states where 

20 the drive frequency and interwell transition rates are much slower than intrawell relaxation 
rates. The signal to noise ratio (SNR) of such a superposition system is given by: 
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SNR = 



n 



+0 



4 l+exp[^o/^2j7i] 



Where Kb is the Boltzmann's constant, T is the temperature, and the sinusoidally 
modulated asymmetry energy e is given by: 

£ = ^0 + xcos^^;^/ 

5 

For the above bistable system, the power spectrum S((o) may be represented by the Fourier 
transform of C(x), containing a roughly Lorentzian broadband noise background and 6- 
fonction peaks at co = 0, the driving frequency cOs, and its harmonics. A measured 
correlation function for the quantum noise is given by C(t) = <n+i (t)n+i (t+ t)>, where each 
1 0 n+i = 0 or 1 , and the <, indicates an average over t over many data points i taken at equal 
intervals, and is given by 

The probability of being in the + quantum state at / after being in the state qo is 
given by 



15 A custom Hamiltonian, couples the above system to an ensemble of harmonic 

oscillators is given by: 



H 
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In the above Hamiltonian, e denotes the asymmetric energy, A is the tunneling 
matrix element, and Ci are Pauli spin matrices, and C,^ is a harmonic oscillator creation 
5 operator with frequencies C0r|» 

An important aspect is to couple the transformed and preconditioned discrete 
microarray output to a mathematical model for a quantum-mechanical dynamical system 
with specific properties. Specific exemplary parameters for use in calculating the 
Hamiltonian are those proposed by A.J. Legett et al., Reviews of Modern Physics, 59, 1, 
10 1 987 and A.O. Caldiera and A.J.Legett Annals of Physics, 149, 374, 1983, The parameters 
are important for an offline simulation of this Spin Boson system on a digital computer. 
The empirical observables are then collected and used to estimate and compute spectral 
properties, which are actually used by the method. 

Experimental or analytical valid parameters for the above system will work with 
15 the technique because the robustness of the method depends only on the bulk and 
qualitative properties shown by this system and the properties of its power spectral 
harmonics. This is an important because the method is substantially immune to nuances 
and specifies of the actual driving mathematical system. 

The information about the effects of the environment is contained in the spectral 
20 density J(o)) 



-18- 



The harmonic amplitudes determine the weights of 5 spikes of an averaged spectral 
power density in an asymptotic state S'^Cco). £ refers to the coupling strength between the 
spin and boson systems, and Po defines the equilibrium state in the absence of driving 
force. Power Spectral Density (PSD) is a well understood concept is signal processing and 
basis for many computational algorithms. For a spectral signature, partitioned into bands of 
interest, power spectral density is formally as the set of measurements of average power in 
each spectral band, normalized by bandwidths. Formal definition can be found in 
"Cyclostationarity in Communications and Signal Processing", Ed. William Gardner, 
1994, IEEE Press, NY, pages 46-47. 

The power amphtudes rj^ in the mth frequency cornponent of asymptotic state 

space are calculated using 



7jQ,^) = 4;r|/^„(Q,£)| 



and the phase shift is given by 




The analytic for the external force is given by 
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pS^,£) = -X^^h (-imQ,r) 
;^-imQ TV 

The parameters Eq are predetermined and are design specific. Typically, values of 
0.001 and 0.0001 are used for y and Eq respectively . In the above expression determine 

5 the weights of the S spikes of the averaged spectral power density. 

For some embodiments, the QEF is designed by matching the power spectral 
density (PSD) amplitude of preconditioned, extracted events of interest to that of the Spin 
Boson system described above so that stochastic and deteiministic time scales match and 
so that the time scales couple back to noise statistics and degree of asymmetry. The 
1 0 technique employs a fully automated iterative conjugate gradient relaxation method for 
spectral niatching between asymmetric base system and the transformed, preconditioned 
extraction core. The determination of the QEF depends on the specifics of bioelectronics 
substrates used for actual analysis. The method is however geiieralizable to all arrayed 
c enibodiments. In addition, the method is highly scalable to array dimensions (as the 

15 offline design trade-space time does not matter to computational complexity). As the 
system is an overdetermined coupled system, convergence criteria and stability of 
relaxation method directly relates to downstream resonance effectiveness. 
Generation of Order Function 

The order function (OF) of ground truth is generated as follows. The order 
20 function (OF) is for ground truth wherein ground truth represents a state wherein a positive 
signal to noise ratio (SNR) is expected for pixel intensities of selected exemplars used for 



-20- 



developing an QEF. 



A classical Order Function (OF) is used with the following form: 



1 

where i = V-I^ and = ^d0h{e)e-^'^''^ and 0 < ^ < 1 , 



The term 



1 

= jd0h(0)e~^'"^^ is referred to as the Diado integral. 



10 In this application, QEF are designed from a canonical dynamical system drawn 

from a family of exponential systems. Daido integral is a specific term used to construct an 
. Order Function, which is then used to construct the QEF. 

The OF is calculated using an order match which implies that the variance of the 
density function for a specific exponential family approaches 0 (within 0.000001 - 0.0001). 

1 5 The free energy for the density function is given by 



where and represent state vectors. represents the random observable in 
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symmetric bilinear form, and X denotes the characteristic function, calculated using 
extractions from exemplars with reference objects of interest known (see D.C. Brody and 
L.P. Hughston, "Geometry of Quantum Statistical Inference", Physics Review Letters, 
77(14), pp. 2851-2855, 1996). 
5 Determination of Entrainment States 

The OF derivation uses results from Daido's theory of multibranch entrainment of 
coupled nonlinear oscillators, wherein a number of different entrained states co-exist For 
exemplars, an a priori measurement of entrainnient is performed. This is done by using an 
approximate Daido Integral() with Zk = PSD maxima at regions where boundary of event 
10 of interest > desired detection threshold (see Physics Review Letters, 77(7), 1406-1411.) 
Pre-multiplier constants are used to ensure that Zk meets the following tests: 

- maxima power spectrum density (PSD) matches; and 

- L2 norm on even harmonics < e where e = NS-MRF Barrier assuming complete 
spatial randomness (CSR) assumptions discussed in U.S. Pat. No. 6,136,541, entitled 

15 "Method and Apparatus for Analyzing Hybridized Biochip Patterns using Resonance 
Interactions Employing Quantum Expressor Functions all of which is hereby fully 
incorporated by reference". Note that NS-MRF refers to the nonstationary Markov random 
field representation. 

The notion of entrainment states is exploited, in part, because the method treats the 
20 spectral system embodied in an extraction core corresponding to the event of interest, as a 
special case of coupled nonlinear oscillators in equilibrium. However, due to device 
imperfections, signal degradation and other limitations, number of entrained binding states 
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(i.e., background and clutter conditions) may coexist. A PSD match is desired where the 
OF is Z-peaked (but single peaked around each event of interest centroid). Absence of a 
single peak implies imperfect or lack of presence of an event of interest. It also defines the 
resonance loci for this method (i.e., where maximum SNR enhancement) is obtained. 
5 Ground Truth based Modulation of OF 

The OF of ground truth is modulated to yield the QEF as follows. Under controlled 
calibration, as stated above, maximal SNR enhancement (optirnal resonance) is achieved 
when OF yields a single peak. It is an important design point for matching PSD of 
coupling Spin Boson system to the synthetic QEF. The specific form of the QEF to be 
1 0 used is the parameters of a generic OF shown above. So the exemplary method exploits 
two connotations of OF: (a) parametric form for the QEF (that is closer to the classical 
forni) and (b) as exponential attractor for a dissipative system. The two OF's are then 
recoupled and convolved with spectrally transformed exemplar used for constructing the 
QEF. 

1 5 The resulting QEF is given by: 

QEF^,,^^^(^) = -|; XL ^ 

Typically, only first three to twenty-five PSD peaks are considered for spectral matching. 

Preferably the QEF is represented digitally using a matrix or array having the same 
number of elements in the spectrally transformed and preconditioned, extraction core to be 
20 analyzed. In the exemplary embodiment this pertains to the preprocessing done to the post- 
hybridization microarray intensity output to get it into a form where it can be convolved 
with the Quantum Expressor Function. 



-23- 



A major limiting restriction in QSR that is avoided by the exemplary method 
pertains to matching the stochastic and deterministic time scales in "domain system" and 
the external coupling asymmetric dynamical system, since this has liniited applicability to 
5 continuous data. By replacing the time scale match requirement with ensemble spatial 
statistics (statistical mechanics of the signal emitters responsible for generating the pattern 
of interest) expressed via the aforementioned generators an entirely new class of analysis 
for discrete systems is enabled. By deriving the condition that signal emitters need to 
satisfy to maximize QEF generator capacity, a properly constructed and conditioned QEF 
10 will be able to extract and enhance an entire class of interest to produce one or more 
expressor fiinctions represented mathematically as 




Preconditioning 



15 



The signal pattern output by platform array detector 102 is represented 
mathematically as: 



I{N,M) 



20 
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Using the input signal pattern, the Preconditioner unit 104 performs the 
mathematical calculations detailed in U.S. Pat. App. Ser, No. 10/430,664 entitled "Method 
and System for Characterizing Microarray Output Data" filed on May 5, 2003 which is 
hereby incorporated by reference, to generate the preconditioned output signal pattern 
5 represented mathematically as: 

10 c. Convolution Coupling 

Active interferometric coupler 106 performs the following mathematical 
calculations using the preconditioned signal pattern: 



20 



r 



is defined as a vector containing the preconditioned components from an event 
15 of interest, 

and y "~ / 



where Q^^ represents is the QEF after / convolutions. 

^(1) -T(0)T:a) 



Thus f f Q 

r -7''^' =r'Q'''Q'" 

where ^ represents the QEF developed in the preceding step and (its dimensipnless 
25 quantity) 
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— (0 

O 

and where ^ represents the i-th perturbation to the QEF, induced by perturbing one 
of its spectral components. 

for k = 1 to 

for j= 1 to 1000 (set to a large counter value) 

perturb the component of QEF as below 



& (k) = (k) + jq siniwj + Q )T 

where 



W ={ } 
•0 if x<0^ 

and 

3 360 ' ^ denotes a small constant; 
Let ^0 ^ the variance computed from the values 



/ 



]\^Clx(^y ) where fpc denotes the preconditioned spectral vector 

corresponding to a known event of interest present in the arrayed pattern 
being analyzed. 

2 

As an example, fpc refers to the spectral components of the positive control 
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The convolution iteration can be expressed as: 




where f^c refers to the spectral components of a canonical negative control, or 
preconditioned footprint of an event of interest known to be absent in the arrayed 
image. 

After each convolution iteration check for monotonicity of 



if yes — > then exit loop to perform global QEF iterations (that means this 

particular k component is impoitant, i.e we are diverging from 
the negative control.) 
If no — > then continue 
endyloop 
end k loop 

Global QEF iterations may be performed if monotonic divergence is detected 
between the preconditioned extraction core being analyzed and the canonical negative 
control, then the same convolution coupling operations are repeated for all the spectral 
harmonics. The global QEF iterations are provided by: 



R 



-kj+2 




>1 AND 



>1 



R 



A 



Set 



-0 

r 





1 to 25 (chosen to be a small count value) 
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f Nc 

where ^1 captures the variance of the components of ^^Qxi^f ^ ) 



and 



C2 =Q+^ 



^Parseval Avg from Pos- Con. PM^ 
Parseval Avg from Neg. Con PM j 



where Parseval Avg. from Pos. Con. PM refers to the parseval number for a canonical 
event of interest known to be present, and Parseval Avg. from Neg. Con. PM refers to the 
parseval number for a canonical event of interest known to be absent. 

and ^ is chosen to be small, 0.0001 . 

Again after each coupler iteration compute the term 

R R R 

Successively compute m m+l ^ m^2 * * * 
After each convolution iteration check if 

^^>1AND^>1 

If the conditions of the above test are met, resonance is concluded and event of interest is 
called present. 
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If the monotonicity test fails, them the preconditioned test pattern is normalized using the 
expressions below. 



pc J Nc for any component 



U) - 70) 



then 



r 



25 



'o-i) ^ 



The above assumes an object with 25 spectral harmonics of interest. 



The detailed equations for the coupler unit and resonance detector unit are given below: 



for 



y = 1 to N 



A7^7 + l = 



-29- 



j *k 



i = k 



A 



2 -^((Pc/J^H^oJ-^^c) 

i = k 



10 



^ J-\J+2 
Check j^J 



>1 



>i 



, 1 2^ 
0 = 

3 360 



where 



B = fiw,) 



W, = (7 



pm,njc 



and 



\ 



Jman PSD 



pm,njc J 
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^ = U + ' {^H^o U +^)^0c )} 

B = 0^r)\sin{wXr) + (l)^,)] 



^1 



{=0.0001} 



W, = a 



PSD 



QEF pc,pm 



V 



max(PSD 



QEF;pc,pm ) J 



10 



W,=CT 



PSD 



QEF,NC,pm 



max(PSD 



QEF,NC,pm) J 
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d. Resonant Marker Identification 

Finally, resonance marker detector 108 performs the following 
mathematical calculations using the convolved signal pattern to identify the events 
of interest within the convolved signal pattern. 



20 
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The resonant iteration is terminated when 



>1 



>1 



5 



or when iteration counter / exceeds preset "iV" (e.g., 10^ iterations) (for digital 



approximation to analog dynamics). 

Preferably, the resonance interaction is performed digitally by applying a matrix 
representative of the resonance equation to a matrix representative of the resonance 

V 

1 0 stimulus in combination with a matrix representative of the preconditioned extraction of 
event of interest. 

The final result of interferometric analysis is the readout of locations wherein 
resonance has occurred and as identified by row and column number (y). 
Resonance Output Interpretation to Identify Gene Expression 
1 5 For example, with regard to gene expression, once resonances peak are observed by 

applying a QEF to specific object extractions, the object addresses (k,!) of those locations 
are mapped into the physical feature or gene table (which contains the complimentary 
RNA or cDNA probe sequences associated with object locations) to thereby identify the 
gene being expressed, if any, present in the sample being analyzed. This is a simple table 
20 look-up resulting in a direct readout of the names of the expressed genes. 
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Tq summarize the open loop interferometric technique of FIGS. 1 and 2, the 
5 preconditioning of arrayed data to the spectral domain is performed independently prior to 
coupling with the expressor function. Hence, the reverberation convergence perfornied by 
the active interferometric coupler to achieve a resonance state is performed independent of 
preconditioning. The reverberations are an open loop process terminated upon detection of 
a predetermined condition (i.e. observation of a resonance marker.) The active 
10 interferometric coupler implements either a destructive or a constructive interference of 
two dynamical systems. 

B. GOAL-DIRECTED INTERFEROMETRIC SYSTEM 

FIG. 3 illustrates exemplary system components for the second exemplary 

15 implementation of the invention. As before, components employed only during a design 
phase are shown on the left and components employed during the actual analysis of an 
input signal are shown on the right. Some of the components of the system of FIG. 3 are 
the same as those shown in FIG. 1 and hence will be not described again in detail. Again, 
during a design phase, an expressor fiinction generation unit 200 generates expressor 

20 functions based upon calibration data and based on the specific characteristics of a 
platform array detector 202 and during the operational phase, the platform array detector 
generates a signal pattern for analysis. However, no separate preconditioner unit is 
provided. Nor is any separate resonant maker detector provided. Rather, a single adaptive 
interferometric coupler 204 is provided that operates to simultaneously precondition the 

25 detected signal pattern to convert the arrayed signal pattem to a spectral domain while 
convolving the signal pattern and the expressor functions to interferometrically enhance of 
portions of the signal pattern and to identify events of interest, if any, within the 
preconditioned signal pattern. Again, with the exception of the platform array detector, 
which is a physical hardware component, all other components illustrated in FIG. 3 may 
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be implemented in software, hardware, firmware or some combination thereof. 

Note that adaptive interferometric coupler . 204 does not merely perform the 
operations previously described with respect components 104, 106 and 108 of FIG. 1. 
Rather, adaptive interferometric coupler 204 performs a different set of operations. 
5 Briefly, with the technique of FIG. 3, couphng intelligence is nierged with preconditioning • 
and the expressor functions are pre-coupled to the preconditioned functions, as will be 
apparent from the detailed example below. Convergent reverberation is performed 
simultaneously with preconditioning, which is now a closed loop process. Hence, with the 
technique of FIG. 3, preconditioning does not merely transform spatial arrayed data to 

10 spectral data but also transforms the data to a pre-determined state of a predetermined 
dynamical system. If that state is reached, then resonance is implicitly assumed to have 
occurred. Thus, completion of the adaptive interferometric coupling procedure implies the 
detection of resonant markers and so separate resonant marker detection unit is not 
required. ' 

15 Method steps associated with the goal-directed interferometric system will now 

briefly be sumnriarized with reference to FIG. 4. At step 206, a signal pattern derived from 
a particular platform array detector is input and, at step 208, expressor functions designed 
to extract spectral invariants for events of interest detectable by the particular platform 
array detector are input Then, at step 210, the input signal pattern is preconditioned to 

20 convert the signal pattern to the spectral domain while convolving the signal with the 
expressor functions so as to emulate an active interferometric enhancement of portions of 
the signal pattern so . as to identify events of interest therein. As with the technique 
described above, a list of specific events of interest need not be input during the 
operational phase. Rather, the technique operates to detect possible events of interest 

25 present in the input pattern based solely on the input pattern and the previously-generated 
expressor functions. 

One of ordinary skill in the art will appreciate that the equations provided above in 
connection with the open loop interferometric system may be accordingly adjusted such 
that the resonant marker unit and the convolution coupling unit may be combined so that 
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the threshold check for resonance is performed during the post-convolution normalization. 
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C. SELF-ORGANIZING INTERFEROMETRIC SYSTEM 

FIG. 5 illustrates exemplary system components for the third exemplary 
implementation of the invention. Unlike the systems thus far described, no separate design 
5 phase is required. Rather, all components operate during the actual analysis of an input 
signal. Again, some components are the saine as those described above and hence will be 
not described again in detail. In use, a platform array detector 302 generates a signal 
pattern for analysis. The signal pattern is preconditioned by a preconditioner unit 304, 
which, as before, operates to convert the signal pattern to a spectral domain in which 

10 spectral harmonics parameterize events of interest to a predetermined dynamical system. 
The preconditioned signal pattern as well as characteristics of the platform array detector, 
calibration data and a set of canonical expressor functions are then all input by an 
expressor function adaptation unit 306. The canonical expressor functions are generalized 
expressor functions that have not been tailored to the particular platform being used. No 

15 separate design phase expressor function generation unit is. provided. Nor is a separate 
resonant maker detector provided. Rather, expressor function adaptation unit 306 operates 
to simultaneously generate preconditioned expressor functions based on the preconditioned 
signal pattern and based on the canonical expressor functions while convolving the 
preconditioned signal pattern and the preconditioned expressor functions to 

20 interferometrically enhance portions of the signal pattern to identify events of interest. As 
before, with the exception of the platform array detector, which is a physical hardware 
component, all other components illustrated in FIG* 5 may be implemented in software, 
hardware, firmware or some combination thereof. 

Note that expressor function adaptation unit 306 does not merely perform the 

25 operations previously described with respect to components 100, 106 and 108 of FIG. 1. 
Rather, with the technique of the self-organizing interferometric system, the- coupling 
intelligence and the reverberation dynamics are pre-built into the expressor functions and 
so the expressor functions are not merely a representation of a spectral marker as with the 
technique of FIGS. 1 and 2 and with the technique of FIGS. 3 and 4. This will be more 
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readily apparent from the detailed example below. Also, although the preconditioning unit 
may be the same as in the technique of FIGS. 1 and 2, the preconditioned array now 
serves as a preconditioner to the expressor function to thereby produce a properly 
preconditioned expressor, which itself corresponds to the detection of resonant marker. 
5 Hence, generating and outputting the properly preconditioned expressor function results in 
detection of the event of interest. As the expressor functions are tailored to the 
preconditioned array, they cannot be used again to process other signal patterns detected by 
the platforrn array detector. Also, as will be described in greater detail below, with the 
technique of the self-organizing interferometric system, the operations of the expressor 

1 0 function adaptation unit operate in an internal closed loop. 

Method steps associated with the self-organizing interferometric system will now 
briefly be summarized with reference to FIG. 6. At step 308, a signal pattern derived from 
a particular platform array detector platform is input and, at step 310, the signal patterns 
are preconditioned to convert the signal pattern to the spectral domain in which spectral 

1 5 harmonics parameterize events of interest to the aforementioned predetermined dynarnical 
system. At step 312, preconditioned expressor functions are then generated based on the 
preconditioned signal pattern and on the characteristics of the platform array detector and 
based on the events of interest. This is performed using a closed loop process employing 
reverberant convergence so as to emulate active interferometric enhancement of portions 

20 of the original signal pattern associated with the events of interest. Hence, generation of 
the preconditioned expressor functions allows immediate detection of the events of interest 
without further processing. As before, no list of specific events of interest rieed be input 
during the operational phase. 

One of ordinary skill in the art will appreciate that the equations provided above in 

25 connection with the open loop interferometric system may be accordingly adjusted such 
that the modulation of the OF to generate a QEF is immediately succeeded by the 
convolution coupling. 
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D. ITERATIVE INTERPEROMETRIC SYSTEM 

FIG. 7 illustrates exemplary system components for the fourth exemplary 
implementation of the invention. As before, components employed only during a design 
phase are shown on the left and components employed during the actual analysis of an 
5 input signal are shown on the right and components already described will be not described 
again in detail. During the design phase, an expressor function generation unit 400 
generates expressor functions based upon calibration data and based on the specific 
characteristics of a platform array detector 402 and, during the operational phase, the 
platforrn array detector generates a signal pattern for analysis. The signal pattern is 

10 preconditioned by a preconditioner unit 404, which, as before, operates to convert the 
signal pattern to a spectral domain in which spectral harmonics parameterize events of 
interest to a predetermined dynamical system. An iterative interferometric coupler 406 
convolves the preconditioned signal pattern and the previously-generated expressor 
functions so as to interferometrically enhance portions of the preconditioned signal pattern 

1 5 associated with events of interest that are present within the preconditioned signal pattern. 
However, unlike the technique of the open loop interferometric system of FIGS. 1 and 2, 
wherein the convolution coupler is a stand-alone unit that operates to generate a final 
convolved signal pattern, with the technique of FIG. 7, the iterative interferometric coupler 
operates under the control of an adaptive controller 408 to sequentially produce an 

20 iteratively convolved signal pattern. As explained more fully below, the adaptive 
controller operates to control the interferometric coupler to iteratively and selectively 
convolve expressor functions to the preconditioned signal pattern until a predetermined 
degree of convergence is achieved so as to allow identification of events of interest within 
the enhanced signal pattern. Again, with the exception of the platform array detector, all 

25 other components illustrated in FIG. 7 may be implemented in software, hardware, 
firmware or some combination thereof. 

Method steps associated with the iterative interferometric system will now briefly 
be summarized with reference to FIG. 8. At step 410, a signal pattern derived from a 
particular platform array detector is input and, at step 412, expressor functions designed to 
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extract spectral invariants for events of interest detectable by the particular platform array 
detector are input Then, at step 414, the input signal pattern is preconditioned to convert 
the signal pattern to the spectral domain spectral domain in which spectral harmonics 
parameterize events of interest to the aforementioned predetermined dynamical system. At 
5 step 416, the preconditioned signal pattern and the expressor functions are then iteratively 
convolved under the control of the controller to achieve a reverberant convergence so as to 
emulate active interferometric enhancement of the signal pattern associated and to allow 
detection of events of interest. As before, no list of specific events of interest need be 
input during the operational phase, 
1 0 One of ordinary skill in the art will appreciate that the equations provided above in 

connection with the open loop interferometric system may be accordingly adjusted such 
that multiple QEFs may be employed to achieve a desired resonance level and the most 
relevant QEF is chosen after each convolution if a desired resonance is not achieved. 

15 E, GENERAL FEATURES 

Exemplary embodiments are described above wherein arrayed data generated by a 
platform array detector is processed to identify events of interest therein. A wide variety of 
types of arrayed data maybe processed. FIG. 9 lists specific examples of types of arrayed 
data that may be processed: The examples shown include spatial 2-D data, spatial 1-D 

20 data, point emitter data, temporal point emitter data, spatio-temporal point emitter data and 
spectral point emitter data. The different types of arrayed data may be separately processed 
or, as shown, combined to yield virtual arrays for processing. In FIG. 9, the symbol \(f 
generally represents a virtual array combiner operation, wherein individual spatial values 
are collected together into a single function or data array. For example, if sixteen 

25 individual point emitters are provided, the data from the sixteen point emitters is combined 
or merged to yield a single 4X4 virtual array (represented, e.g., as f(x,y)). The symbol i 
generally represents a temporal accumulation operation, wherein individual values are 
accumulated over time to yield a single time-varying function or data array. For example, 
if a single point emitters is provided, the point emitter may be sampled periodically to 
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accumulate values to populate a time-varying function (represented, e.g., as f(t)). In some 
cases, both virtual array and accumulation operations are performed to generate spatio- 
temporal data (represented, e.g., as f(x,y,t)). 
More specifically, FIG. 9 illustrates: 
5 (a) spatial data from a 2-D spatial array 500 represented in the form of a 2-D 

function h(x,y); 

(b) spatial data from a 1-D spatial array 502 represented the form of a 1 -D function 

k(x); 

(c) spatial data from a group of point emitters 504 (with each point emitten, for i == 
10 1 ... n, providing a single non-time- varying value) assembled into a single 2-D fonction 

h(x,y) via a virtual array combiner 506; 

(d) temporal data from a single point emitter 508 accumulated over time via a 
temporal accumulator 5 1 0 to yield a single l-D function g(t); 

(e) spatio-temporal data from a group of individual point emitter arrays 512 (with 
15 each point emitten, for i = 1 ... n, providing a 2-D anay non-time varying values) each 

separately accumulated over time via temporal accumulators 514 to yield individual 
functions p(n,y,t), which are in turn assembled into a 2-D spatio-temporal fiinction via a 
virtual array combiner 516; and 

(f) spectral data from a group of individual spectral point emitters 518 (with each 
20 spectral point emitteri, for i = 1 ... n, providing an output function f(X) varying with 

wavelength) assembled into a 3-D virtual spectral array via virtual array combiner 520, 
with values then accumulated over time via temporal accumulators 522 to yield a 4-D 
function g(x,y,A,,t). 

Additionally, the various types of arrayed data ((a) - (f)) may be further combined 
25 into a larger virtual array (g) by virtual array combiner 524. Also, any of the types of 
arrayed data ((a) - (g)) may be further filtered via linear or non-linear filters 526, to yield 
filtered array fiinctions (h). The final collection of data is represented in the figure as 
arrayed data 528. As can be appreciated, a wide variety of arrayed data may be input, 
combined and manipulated to yield individual multi-dimensional input functions for 
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processing by the above-described systems. 

Also, a wide variety of preconditioner functions may»be employed for use in 
preconditioning the arrayed data and a wide variety of expressor functions may be 
generated for use in processing the preconditioned arrays. FIG* 10 lists specific examples 
5 of types of expressor functions and types of preconditioning fiinctions. The types of 
expressor functions ^ are identified in terms of the type of noise associated with the 
expressor function. Examples of expressor functions include quantuni expressor functions 
(QEFs), classical expressor functions (CEFs), classical statistical noise expressor functions 
(SCNEFs), pseudorandom noise expressor functions (PRNEFs), systemic bias expressor 

10 functions and expressor functions based on polycyclostationary noise (PCS), 
cyclostationary noise (CS), stationary noise (SN), non-stationary noise (NS). Examples of 
preconditioner furictions include 1-D Fourier functions, 2-D Fourier functions, N-D 
Fourier functions, a time division multiplexing (TDM) functions, wavelength division 
multiplexing (WDM) functions, frequency division multiplexing (FD) functions, radial 

15 basis functions, wavelet kernel functions, fractal functions and soliton functions. In 
general, any combination of expressor function type and preconditioner function type 
illustrated in the figure may be employed. 

FIG. 11 illustrates an exemplary system 700 incorporating a biological platform 
array 702 that provides arrayed data to a software -based active interferometric analysis 

20 system 704. The platform array includes a solid support 706 upon which various probes 
708 are attached via attachment chemistry portion 710. In the example of FIG. 11, the 
probes include one or more of cDNA, mRNA or oligonucleotide probes or proteins, 
peptides or oligopeptides. Following application of a biological sample, such as blood or 
urine, to the probes, a detector 712 detects a resulting pattern, which is then forwarded to 

25 the analysis system in the form of arrayed data. The analysis system exploits one of the 
general techniques described above to identify events of interest present in the biological 
sample. In this example, the events of interest are typically output in the form of one or 
more "gene calls". 

A wide variety of different platform array detectors may be used to generate the 
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arrayed data. Examples of platforrns detectors are listed in FIG. 12. The examples include 
optical platforms, biomolecular platforms, ionic platfonns, biomechanical platforms, 
optoelectronic platforms, radio frequency platforms, and other electronic microdevices. 
Exemplary biomolecular spatial array platforms include: hybridized spotted cDNA 
5 microarrays, synthesized oligonucleotide arrays, spotted oligonucleotide arrays, peptide 
nucleotide assays, single nucleotide polymorphism (SNP) arrays, carbohydrate arrays, 
glycoprotein arrays,, protein arrays, proteomic arrays, tissue arrays, antibody arrays, 
antigen arrays, bioassays, sequencing microarrays, sequencing by hybridization (SBH) 
microarrays, siRNA duplexes, RNAi arrays glass-based arrays, nylon membrane arrays, 

10 thin film arrays, polymer-substrate arrays, capillary electrophoresis arrays, genospectral 
arrays, electronic arrays, bead arrays, quantum dot arrays, glycan arrays, spotted wells, and 
spotted well plates. Typically, any of the foregoing platforms may be implemented as a 
microarray. Components of a typical microarray are shown in FIG. 13, which include a 
substrate, a probe portion, a label portion, a bioassay portion, and a detector. Exemplary 

1 5 substrates include glass, nylon, thin film, and polymer. Exemplary probes CDNA probes, 
spotted oligonucleotide probes and synthetic oligonucleotide probes. Exemplary labels 1- 
dye labels and 2-dye labels. Exemplary bioassay portions include prehybridization buffers, 
hybridization buffers and wash buffers. Exemplary detectors include laser scanners, 
confocal microscopy detectors, charge coupled devices (CCDs) and electronic readouts, 

20 The substrate has a certain layout represented by feature size, various controls (positive, 
negative dye and alignment), and probe density. 

FIG. 14 provides a summary of an exemplary technique of the invention. Briefly, 
the figure illustrates method for actively analyzing a signal pattern representative of 
arrayed data to identify events of interest therein. At step 700, a signal pattern 

25 representative of arrayed data is input and, at step 702, resonance patterns are generated 
based on interference between synthetic noise and the signal pattern, in accordance with 
techniques already described. Thereafter, at step 704, resonances are detected within the 
resonance patterns associated with events of interest. The synthetic noise used at step 702 
is preferably in the form of a quantum expressor function, a classical expressor fimction. 
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classical statistical noise, pseudorandom noise, a systemic bias or some combination 
thereof. The arrayed data may be, as already explained, in the form of a spatial 2-D array, 
a spatial 1-D array, an N-D array, a temporal point emitter array, a spatio-temporal point 

emitter array, spectral point emitter array or a virtual array constructed by combining 

J 

5 spatial separate point emitters. 

The description of the exemplary embodiments is provided to enable any person 
skilled in the art to make or use the present invention. Various modifications to these 
embodiments will be readily apparent to those skilled in the art and the generic principles 
defined herein may be appUed to other embodiments without the use of the inventive 
10 faculty. Thus, the invention is not intended to be limited to the embodiments shown herein 
but is to be accorded the widest scope consistent with the principles and novel features 
disclosed herein. 



